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Abstract 

We explore spherically symmetric solutions to the Einstein-Klein-Gordon equations, the 
defining equations of wave dark matter, where the scalar field is of the form /(t, r) — e*"*F(r) 
for some constant w G M and complex- valued function F(r). We show that the corresponding 
metric is static if and only if F{r) — ft,(r)e*" for some constant a £ M and real-valued 
function h{r). We describe the behavior of the resulting solutions, which are called spherically 
symmetric static states. We also describe how, in the low field limit, the parameters defining 
these static states are related and show that these relationships imply important properties 
of the static states. 



1 Introduction 

The study of dark matter in the universe has been an exciting field of research in physics for 
decades and many notions and ideas have been presented to describe this exotic matter. One 
idea is using a scalar field to describe dark matter. This idea has been considered for at least 
twenty years, having been studied under the name of boson stars or scalar field dark matter 
[l}|3|[6|[7| [To|[T2|jl4|[l6|[l8[|2l] . In addition, Bray f4j has presented geometrical motivations for 
describing dark matter in this way. 

The Einstein-Klein-Gordon equations (see equation ([T])) are the defining equations of this 
concept of dark matter. Due to the fact that the Klein-Gordon equation (see equation ( |lb| )) is 
a wave-like equation, we refer to this notion of dark matter as wave dark matter. 

Describing the Einstein-Klein-Gordon equations in spherical symmetry is an excellent place 
to start studying wave dark matter for two reasons. First, it immensely simplifies the equations 
involved making them easier to solve. Secondly, as with many theories, there are some very 
interesting objects to be described that are approximately spherically symmetric. In the case 
of wave dark matter, dwarf spheroidal galaxies are approximately spherically symmetric and 



almost entirely dark matter 15 ,22 implying that their kinematics are controlled by their dark 
matter component. Thus determining what the wave dark matter model predicts in spherical 
symmetry would be important in showing the level of compatibility of wave dark matter with 
dwarf spheroidal galaxies. 

The purpose of this work is to present a few results in this regard. In particular, we describe 
an important class of spherically symmetric solutions of the Einstein-Klein-Gordon equations 
called static states and discuss some properties of these solutions. We should also reference here 



'Mathematics Department, Duke University, Box 90320, Durham, NC 27708, USA, alrparry@math.dulce.edu 



1 



A. Parry 



Spherically Symmetric Static States of Wave Dark Matter 



others' results that, with the exception of the ground state, these static states are unstable under 
perturbations [l}[2||5H8}[ll|[l9], which poses a problem in using the static states as a physical 
model. However, in physical situations dark matter is always coupled with regular matter, which 
may have stabilizing effects on the dark matter. Thus understanding the stability of these static 
states in the presence of regular matter is an important open problem. It is also possible that 
dark matter may not exist as a single static state at all and as such, finding other kinds of stable 
configurations of wave dark matter is another important open problem. These open problems 
are not addressed here but are left for future work. 



2 The Spherically Symmetric Einstein-Klein-Gordon Equations 

In this section, we present the form the Einstein-Klein- Gordon equations take in spherical sym- 
metry. This has been surveyed in detail by the author in [17] and we refer the reader to that 
paper for the derivation of the equations presented here. We also note that several other ref- 
erences including, but certainly not limited to, [l}{3|[5||8|[9|[Tl|[l2}[l4}[l6|[l^ have written the 
Einstein-Klein-Gordon equations in spherical symmetry using either the metric presented here 
or another form of a general spherically symmetric metric. 

To begin, let (A^, g) be an asymptotically Schwarzschild spacetime whose metric has signature 
( — h ++). Let / : — 7- C be a smooth complex scalar field on the manifold N. Then the 
Einstein-Klein-Gordon equations are 



df (g) df + df (E> df \df\ 



G = l^ --' - + |/|2J gj (la) 

Ogf = TV (lb) 

where G is Einstein curvature tensor and Og is the Laplacian with respect to the metric g. The 
parameter T is a fundamental constant of this system and its value is important in using this 
system as a model of dark matter in the universe. The other parameter /xq is not fundamental 
to the system and can be completely absorbed into / if desired. 

In [4], wave dark matter is modeled with a real-valued scalar field. The system ([T]) is the 
complex version of the Einstein-Klein-Gordon equations described in [4] and the equations here 
reduce to the ones in [4j by simply requiring that / is real-valued, i.e. Im(/) = 0. We use the 
complex version here because, as we will see, it possesses static solutions. 

We next impose spherical symmetry. We use the Newtonian-compatible metric surveyed in 
[it], namely. 



for real- valued functions V and M. As surveyed in [17j , the function M can be interpreted as 
the mass of the system, being the Hawking mass of a metric sphere and the flat integral of the 
energy density of the system. By requiring the low field limit, that is, M(t, r) <C r for all t and 
r, we also get that V is approximately the gravitational potential of the system. In this metric. 
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the Einstein-Klein- Gordon equations produce the fohowing set of PDEs [17^ 
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(3c) 
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The function p is defined by equation (3c). There are two other equations produced by the 
Einstein-Klein-Gordon equations which are automatically satisfied by solving the system Q 



later. 



17 . We list them here for completeness and because the first of the two equations will be used 
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We note next the important behavior of these functions at the origin and as r — )• oo. 

We use the previously mentioned fact that M is the flat integral of the energy density to shed 
light on the initial behavior of the function M near the origin, r = 0. Using the construction in 



17 , let ft be the unit vector in the dt direction and let T denote a stress energy tensor, then the 



energy density of this system is defined as fi{t,r) = T{vt, ut), and the Einstein equation yields 



M{t,r) 



I ^{t,s)dV= I ATTs'^fi{t,s)ds 

J Br Jo 



(5) 



Since the energy density is spherically symmetric and smooth, if the energy density is nonzero 
at the central value, r = 0, at any time t, then the energy density must be an even function of 
r. Moreover, since A*(t,T) is also smooth at the origin, fir{t,0) = for all t. Thus for small r, ^ 
is approximately constant and nonnegative. Then the above integral yields for small r that 



/•r pT 

M{t,r)= / 'i7rs'^fj.{t,s)dr ^ / 47rsV(^, 0) 
Jo Jo 



dr = ^^^.3 



(6) 



Thus the initial behavior of M near r = is that of a cubic power function. In particular, this 
implies that for all t 



M(t,0) = 0, M^(t,0) = 0, 



and 



Mrr{t,0) = 



(7) 
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This is consistent with the fact that the metric functions, e^^ and ( 1 



2M\ 



r j 



, are also smooth 



spherically symmetric functions that are nonzero at r = and hence even functions of r. This 
implies that, since r is an odd function, M(t, r) must also act like an odd function near r = 0. 
Similarly, r) must be an even function of r, implying that 



K(t,0) = 

for all t. Equation Q implies via L'Hopital's rule that 



lim — 
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(8) 



(9) 



Next, since / and p are spherically symmetric and allowed to be nonzero at r = 0, we have that 
both / and p are even functions in r as well, which, for regularity at r = 0, implies that 



/r(i,0) = 



and 



p,(t,0) = 0. 



(10) 



for all t. 

Next we consider the behavior of the functions at the outer boundary. Since N is asymp- 
totically Schwarzschild, there exist constants, m > 0, called the total mass of the system, and 
K > 0, and a Schwarzschild metric gs given by 
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such that g approaches 55 as r — t- oo. This yields the following asymptotic boundary conditions. 
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(12) 
(13) 



The first boundary condition (12) implies by equation (3a) that — t- as r — )• 00 and hence 
M approaches a constant value as r — )• 00. Given equations ([2]), (11), and the second boundary 
condition (13), this constant will be the parameter m in (11). 

Now that we have described the Einstein-Klein-Gordon equations in spherical symmetry, we 
are ready to discuss the class of spherically symmetric solutions to the Einstein-Klein-Gordon 
equations that yield static metrics. 



3 Spherically Symmetric Static States 

In the remainder of this paper, we will consider solutions of the spherically symmetric Einstein- 
Klein-Gordon equations where the scalar field / is of the form 

fit,r)=e'^'Fir) (14) 

where cj G M and F is complex-valued. As t increases, / rotates the values of the function F{r) 
through the complex plane with angular frequency uj without changing their absolute value. 
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Thus, without loss of generality, we will assume that w > 0, since if u < 0, F wih simply rotate 
in the opposite direction. With / of this form, we have that 

(15) 
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With a scalar field of this form, the system ([s]) implies that the metric is static under certain 
conditions on the function F. We collect this result in the following proposition. 

Proposition 3.1. Let {N,g) be a spherically symmetric asymptotically Schwarzschild spacetime 
that satisfies the Einstein- Klein- Gordon equations for a scalar field of the form in (14)- Then 
(N, g) is static if and only if F[r) = h{r)e^"' for h real-valued and a € M constant. 

Proof. By definition, a spacetime metric is static if there exists a timelike Killing vector field 
and a spacelike hypersurface that is orthogonal to the Killing vector field. For the metric in 
equation ([2]), if the metric components, V and M, do not depend on t, then dt is one such 
Killing vector field and it is already orthogonal to the t = constant spacelike hyper surfaces. 
If a spherically symmetric metric is static, then we can choose the t coordinate to be in the 
direction of the timelike Killing vector field, making dt the timelike Killing vector field in these 
coordinates, and choose the polar-areal coordinates on its orthogonal spacelike hyper surf aces, 
yielding a metric of the form in Q. Then since the metric cannot change in the direction of dt, 
the metric components must be t-independent. It remains then to show that under and only 
under the given conditions on F, the metric components are t- independent. 



Note that by equation (4a), Mt = if and only if Re{frp) = 0. Using equation (15) and 



(|16|), we obtain 

Re.(Ln) = Re ( e'^* F' (r)(^-^ M - ^ , 

r / 



Reifrp) = Re ( e'^^F'{r)e-^ ( 1 - — ^ ^ (-iuje-''^^F{r)' 



e"^ ( 1 - — j Rei-iujF'{r)F{r)) (17) 



Thus Re{frp) = if and only if Re{—iujF'{r)F{r)) = 0, which is true if and only if F'{r)F{r) is 
real-valued. Now any complex-valued function F{r) can be written as 

F{r) = /i(r)e*"(^) (18) 

for real-valued functions h and a. If we write F{r) this way, then F'{r) is as follows. 

F'{r) = /i'(r)e*'*(^') + ih{r)a' {r)e'^^''^ = e'^'-''^ {h'(r) + ih{r)a'{r)) . (19) 

Then we have that 

F'{r)F{?) = h'{r)h{r) + ih{rfa{r). (20) 
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Since h and a are both real-valued, we see that F'(r)F{r) is real if and only if either h{r) = 
or a'(r) = 0. If h{r) = 0, then F{r) = h{r)e^°' with a G M constant is trivially true. If a'{r) = 0, 
then a(r) is constant and the result still holds. Thus Mt = if and only if F{r) = /i(r)e*" with 
h real- valued and a E M constant. 

It suffices from here to show that Mt = if and only if g is t-independent. Obviously, g 
being t-independent implies Mt = 0. On the other hand, assume that Mt = 0. First note that, 
since Mt = 0, Mrt = as well. Moreover, since |/|^ = |-F|^ and Ifrf' = \F'\'^ and F has zero 



t-derivative, then |/p and \fr\'^ both have zero t-derivatives. Differentiating (3a) with respect 
to t and using the fact that l/l'^, |/r|^, and M all have zero t-derivatives, we have that 



Mrt = 47rr^//o 



= 47rr>o 1 
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(21) 



Thus since, M, T, and //q are all nonzero, dt[\p 



0. 



Next, since (3a) and (3b) completely det ermine the Einstein equation, the function V is 
determined at every value of t by solving (3b) at that value of t. Since |/|^, l/rp, |pp, and M 



never change with t, Vr never changes with t and hence the solution, V , of (3b), by uniqueness 
of the solution to an ODE, will never change with t so long as V{t^ 0) = constant. Since V{t, 0) 
is determined to be the value that makes V satisfy (13) and since M never changes with t, 
then the fact that Vr is t-independent forces V{t, 0) = constant. Thus = and the metric is 
t-independent. □ 



Since the value of a simply adjusts the "starting position" of the values of / before they 
rotate, we will, without loss of generality, set a = 0, which is the same assumption that F{r) 
be real-valued. This amounts to simply choosing the hypersurface that we denote as t = 0. 



To summarize, we restrict our attention to static states of the form in (14) with w > and F 
real- valued. 



3.1 ODEs for Static States 



In this section, we input our ansatz (14) with the requirement that uj > {) and F is real- 
valued into the Einstein-Klein-Gordon equations (3a)-(3d). Since the metric is t-independent, 
we summarize that 
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Additionally, if we differentiate (16) with respect to t and (15) with respect to r, we obtain 
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Equation (3c) turns into (16) becoming redundant. The last equation, (3d), becomes 
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which yields 
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To make the system first order, we will introduce a new function H(r) = F'{r). Then (28), (29), 



and (31) become the system 
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3.2 Boundary Conditions 



We will solve the system (32) numerically, but in order to do so, we need to express how we will 
deal with our boundary conditions numerically. Ideally, we would like to model the system in an 
infinite spacetime, but since we are computing these solutions numerically, we must introduce an 
artificial right hand boundary, say at r — '^max'i to which we restrict our domain. To simulate the 
fact that the spacetime is asymptotically Schwarzschild, which we detailed in (12) and (13), we 
will choose r^ax sufficiently large and impose the condition that the spacetime be approximately 
Schwarzschild at the boundary. That is, we will impose (12) and (13) at the boundary r = r^ax- 



In this case, equation (13) at r = r^ax becomes 
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For (12) at r 



we require / to approximately solve the Klein-Gordon equation (lb) 



in the Schwarzschild metric (11). Computing the Laplacian in the Schwarzschild metric, this 
equation becomes 
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This ODE is routinely solved and has the general solution 
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for some constants Ci, C2 G M. However, we also require that F — )• as r — )• 00 so that / — )• 
as well. Thus Ci = and we relabel C2 as simply C. That is, at T — f^max 1 ^6 r6C[Uir6 



F = ^e-^^^. (37) 
r 

We have no way of directly determining the correct value of C associated with a given static 
solution. However, if we differentiate the above with respect to r, we obtain 




(38) 

which does not depend on C. Then the condition that at r = Vmax, f approximately satisfies 
the Klein-Gordon equation with the Schwarzschild background metric reduces to requiring that 



F'irmax) + F{r^ax) ~ 0. (39) 



This condition imposes that F is decaying appropriately to 0. It also puts a restriction on the 
possible values of uj. Since the left hand side of the above equation must be real, we have that 
— ^ > 0, or equivalently, 



UJ 



- < T (40) 

K 



That is, uj/ K G [0, T]. In our numerical calculations, we normalize this quantity and keep track 
of — — E [0, 1] instead. 

Next we consider the boundary at the origin r = 0. We have already noted above that 
M(0) = 0. Moreover, since /r(t,0) = for all t, H{Q) = F'(0) = e-*'^*/,.(t, 0) = 0. Note 



that by the system ([32|, iJ(0) = imphes that if F(0) = 0, then F{r) = for ah r, since H'{r) 
and F'{r) will never change. We are interested in non-trivial solutions to these equations so we 
require that F{Q) ^ 0. However, for any constants c and if c/ is a solution to the Einstein- 
Klein-Gordon equations, ([T]), with /io = Z^*, then / is a solution to the Einstein-Klein-Gordon 
equations with /^o = c'^/i*. Thus, without loss of generality, we will set F{Q) = 1 absorbing any 
excess factors into /^o- 

At this point, we have four remaining parameters to choose, namely, 1^(0), T, and /io- 



By requiring (33) and (39), two of these values are determined, leaving two remaining degrees 
of freedom. The parameter T is an as yet unknown fundamental constant, making it important 
in our computations to be able to freely set T so that we can test different values and see how 
they match the data. The parameter /io controls the magnitude of the energy density and so 
seems a natural choice as another parameter to freely choose. However, this choice is not the 
only choice that could be made. For example, one could instead freely choose T and oj and use 
the boundary conditions to pin down F(0) and no with equivalent results. 

When freely choosing T and /io, to compute the other two parameters, F(0) and lj, we solve 
a shooting problem to satisfy the desired boundary conditions. We illustrate next in detail how 
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we performed this shooting procedure. However, first, we summarize the information about the 
boundary conditions. 

For some choice of T and we require at r = 0, 



F(0) = 1, H{0) = 0, 

and choose ui and Vq to satisfy 



M(0) = 0, 



V{0) = Vo, 



(41) 

(42) 
(43) 

For simphcity in our calculations, we set k = 1. Then we use the standard forward Euler method 




to solve the system (32) with these boundary conditions. 



To understand the procedure we used to solve the shooting problems mentioned above, we 



first comment about what equation (31) tells us about the behavior of F in our system. We 



are solving these equations in the low field limit where the metric ^ is close to Minkowski. 
That is, both V and M are approximately zero. Recall that equation (31) came immediately 
from equation (30), which results from substituting (14) into (3d). Letting V = M = makes 



equation ( 3d ) and the first line of equation ( [SO] ) reduce to 



ArF 



UJ' 



(44) 



where A,, is the Laplacian in 
above equation is 



in spherical coordinates. The one dimensional analogue to the 



kF. 



(45) 



The solution of this differential equation depends on the sign of k. If k > 0, then the solutions 
either exponentially grow or decay. If A; < 0, the solutions exhibit oscillatory behavior. In 



equation (31 ), the analogous coefficient that will control whether the solutions exhibit oscillatory 



or exponential behavior is the following, which we will denote as A(r). 



A(r) 



^2 



(46) 



The sign of A(r) depends on r. While A(r) < 0, the solution of (31) will exhibit oscillatory 
behavior. On the other hand, while A(r) > 0, the solution of (31) will exhibit exponential 
growth or decay. 

Since T and w are constants, the sign of A(r) is completely controlled by the exp(2y) term. 
If V is strictly increasing and starts low enough, then for small r, exp(2y) will be small yielding 
that A(r) < and the initial part of the solution will oscillate. Then as r increases, exp(2V) 
will eventually be large enough that the term dominates A(r) making A(r) > and the end 



behavior will be exponential growth or decay. Imposing boundary condition (42) will ensure 
decaying end behavior instead of growth. 

The two parameters that control where this change from oscillation to exponential decay 
occurs are the initial value of V and the parameter u. Larger values of uj and lower values of 
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V{0) will increase the length of the oscillating region and push the point where it changes to 
exponential decay out to larger radii. In fact, given a fixed T and ^uq, for each value of ^(0), 
there is a discrete infinite set of uj values for which each uj in the set corresponds to a solution 
F that has a given number of zeros (caused by a lengthening of the oscillating region) and the 
appropriate end behavior. 

Thus we perform our shooting problem as follows. To find a solution with say n zeros, first, 
fix a choice of T and hq. Then pick a value of V{0) = Vq. Since we have set k = 1, cj < T in 
order to be able to satisfy equation ( 42 ) . Then we systematically pick different values of u in 
the interval [0, T] until we obtain the appropriate number of zeros and satisfy (42). Finally, we 
use a Newton's method approach to find the value of Vq whose corresponding solution with n 
zeros yields a potential function satisfying (43). 

These solutions, now characterized by the number of zeros the resulting scalar field has, are 
referred to as spherically symmetric state states of wave dark matter. A static state with no 
zeros is called a ground state; with one zero, it is called a first excited state; with two, it is 
called a second excited state, and so forth. There is a considerable amount known about static 
states as they have been studied for decades, see [H[2|[5|[8 11 19] for just a few examples. In the 
remaining sections, we will present some useful results about these static states. 



3.3 Plots of Static States 

In Figures [Tito |4j we have plotted examples of the scalar field F{r) (see Figure [l]), Mass M(r) 
(see Figure [2]), energy density /x(r) (see Figure [s]), and gravitational potential V{r) (see Figure 
|4| for a generic ground state and first through third excited states. We make here the following 
three observations. First, in these plots, V is strictly increasing, as expected, and M approaches 
a constant value as r — )■ 00, also as we expected. Second, for each zero of the function F, there 
is a zero of the energy density plot //, a ripple in the mass profile M, and slight but rapid change 
in the slope of the potential V. And finally, the energy density appears to be approximately 
proportional to \F\'^ = 



4 Families of Static States 

We explained above that four parameters control what static solution is generated by the equa- 
tions. However, since we require two conditions on the boundary, choosing only two of these 
parameters will define a static state. As stated before, we choose to define T and fiQ and solve the 
shooting problem for the parameters Vq and u. For each n, this defines a function, 5„ : — )• M^, 
which maps the pair (T,^o) to the pair (Vo,^^^) such that the choices T, ijo,Vq,uj yields an n*^ 
excited state (n = 0, of course, referring to the ground state). 

A natural question to ask here would be "Is there an expression for 5„ for each n?" The 
answer to this question is yes, at least in the low field limit. In fact, we have also found 
expressions for the total mass m of the system and the half-mass radius rh, that is, the radius 
for which M{rh) = m/2. These expressions were found by numerically computing the states for 
several different values of T and /xq, all of which yield a state in the low-field limit. We analyzed 
the resulting values and found that certain log plots between the values were linear. We have 
collected in Figure [5] some of these plots for the ground state that led to this conclusion. 
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^±. 



Plot of F(r) for a Ground State piot of F(r) for a 1^' Excited State 

1 



-1 




Plot of F(r) for a 2"^ Excited State 



Plot of F(r) for a 3^^ Excited State 





Figure 1: Plots of static state scalar fields (specifically the function F{r) in (14)) in the ground 
state and first, second, and third excited states. Note the number of nodes(zeros) of each 
function. 
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Figure 2: Mass profiles for a static ground state and first, second, and third excited states of 
wave dark matter. 
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Energy Density Profile for a Ground State 



a 
Q 

c 
UJ 



Energy Density Profile for a 2"*^ Excited State 




Energy Density Profile for a 1*' Excited State 




r 



Energy iDensity Profile for a 3'^'^ Excited State 




Figure 3: Energy density profiles for a static ground state and first, second, and third excited 
states of wave dark matter. 
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Potential for a Ground State Potential for a f ' Excited State 




Figure 4: Plots of the potential function, V, for a static ground state and first, second, and third 
excited states of wave dark matter. 
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Figure 5: Left: Log plot of the parameters uj and /j-q for a ground state and constant value 
of T = 5. The slope of this plot is almost exactly 1/2. We get the exact same slope for 
other values of T, thus a;(T,/xo) = Te'~^-f'''=9''="'=«(^)v^. Right: Log plot of the parameters 
C frequency i"^) and T for a ground state. The slope of this plot is almost exactly —1. Thus 
C frequencyC^) = C frequency , where C frequency IS a Constant. Similar plots exist for any n*^ 
excited state. 



These log plots yielded the following expressions, which we emphasize are only expected to 
hold in the low field limit. Let w"', V^, m", rJJ be respectively the values of u, V{0), m, and rh 
for an n*'^-excited state corresponding to a choice of T and /io- Then we have that 



tj"(T, /io) T exp (^C]^equency^^^ (47a) 



y^o'i'^^ f^o) ^ C potential (47b) 

m"(T,/.o)«C;;_T-3/Vy' (47c) 
rUr, M) « Cad»sT-^/Vo (47d) 

for some constants C^^^g^ency C^otentiai^ C^ass and Cl^^aius which depend only on n. We have 
computed these constants for the ground through fifth excited states as well as for the tenth 
and twentieth excited states and have collected their values in Table [T} Note also that equations 



(47a) and (47b) constitute the 5„ function mentioned above. 



In the wave dark matter model, T is a fundamental constant that should be the same 



throughout the universe. In the case of constant T, for each n, the equations in (47) define 
one-parameter families of the n**^ excited states, the parameter being the scaling constant /io- 
For constant T then, we have that these static states only differ in size where, by equations 



(47c) and (47d), a larger fiQ corresponds to a more massive and more dense wave dark matter 
halo. 
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n 


fvcQucncy 


potsntittl 


mass 


Tddius 





-3.4638 ±0.010 


-6.7278 ±0.003 


4.567 ±0.05 


0.8462 ± 0.004 


1 


-3.2422 ±0.012 


-7.5411 ±0.007 


10.22 ±0.10 


2.2894 ± 0.009 


2 


-3.1566 ±0.014 


-7.9315 ±0.009 


15.81 ±0.16 


3.8253 ±0.014 


3 


-3.1062 ±0.015 


-8.1823 ±0.010 


21.37 ±0.22 


5.3994 ±0.018 


4 


-3.0714 ±0.015 


-8.3653 ±0.010 


26.91 ±0.27 


6.9860 ± 0.022 


5 


-3.0452 ±0.016 


-8.5086 ±0.011 


32.42 ± 0.33 


8.5606 ±0.026 


10 


-3.0076 ± 0.052 


-9.0018 ±0.037 


60.32 ± 1.18 


15.1357 ±0.039 


20 


-2.9949 ± 0.077 


-9.5061 ±0.074 


116.62 ±2.57 


29.6822 ±0.107 



Table 1: Values of the constants in the system (47) for the ground through fifth excited states as 



well as the tenth and twentieth excited states. We have given them error ranges which encompass 
the interval we observed in our experiments, but it is possible that values outside our ranges 
here could be observed, though we don't expect them to be so by much if the discretization of 
r used in solving the ODEs is sufficiently fine. Note also that our values have less precision as 
we increase n. This is because as n increases, it becomes more difficult to compute the states 
with as much precision. 



4.1 Scalings of Static States 

Certain scalings exist for various quantities if we scale time and length in any coordinate system. 
In particular, let us scale the time coordinate, t, and the standard spatial coordinates, x*, so 
that 

Time: t = fit Distance: = ax* (48) 

for some positive constants /?, q G M. Then velocities in the (t, x) system will scale to velocities 
in the (t, x) system as follows 

dx* a dx* a i , , , 

= = -— = -v"- (49) 
dt P dt (3 ^ ' 

Similarly, other quantities scale as follows (velocity is included again for completeness). 

_ ol _ oi 

Velocity: v = —v Mass: m = -^m 

1 - a2 

Energy Density: U = Gravitational Potential: V = —kV 

Frequency: A = -^A (50) 

These scalings are routine to derive and follow directly from how the units on each of these 
quantities scale with the scaling for mass being that which is required to keep the universal 
gravitational constant the same from one scaled coordinate system to the other. 



From the system (47), the scalings in (48) and (50), and given that the constants C" are 



dimensionless, we infer how the parameters and T would scale under these coordinate scalings 
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in order to keep (47) invariant. To do this, let ci and C2 be such that under the coordinate scalings 
in pSl), 



and 



T = C2T. 



(51) 



Since m and are a mass and spatial quantity respectively, equations ( 47c ) and ( 47d ) yield 



m 



a 
J' 



T-3/2„V4 



p2 ^rnas 
3 3/2 

r, I mast 



3/2^-3/2 -1/4 1/4 



Y-3/2-1/4 



and 



(52) 



arf^ 

^-1/2-1/4 



- or" 1/2^-1/2 1/4--1/4 

— aCg C]^ ^radius PO 

Then requiring these equations to be invariant under coordinate scalings is equivalent to 



33/2 



a c. 



1 



and 



1/2^1/4 
-I 



Solving these two equations simultaneously for ci and C2 yields 



(53) 



(54) 



ci 



1 



and 



C2 



/3 



Thus equations (47c) and (47d) imply that and T scale as follows 

1 



Stress Energy Tensor Constant: /xq 



/3 



Upsilon: T 



/3 



T. 



(55) 



(56) 



We observe here that ^0 scales as energy density, which is expected given that as we said before 
it controls the magnitude of the energy density function defined by the stress energy tensor. 

Since T is a fundamental constant in this system, it is appropriate to ask under what kinds 
of coordinate scalings of the form in (48) is T invariant, that is, T = T. In light of (56), the 
answer to this question is readily apparent, those scalings where /3 = a'^. Under this type of 
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scaling, (48), (50), and (56) become 

Time: t = 

Velocity: v = 
Energy Density: fl = 
Frequency: A = 
Stress Energy Tensor Constant: /xq 



ah 



a 



1 



Distance: x = 
Mass: m 

Gravitational Potential: V ■ 



ax 

1 

a 
1 

' a2 



-m 



V 



a'* 



Upsilon: T = T. (57) 



The scalings in (57) are also consistent with keeping the remaining relations (47a) and ( |47b| ) 
invariant under coordinate scalings. Showing this for (47b) is straightforward and follows from 
(57) and the fact that Vq has gravitational potential units. 



Vr 







2 ^0 



^potential 



T 



(^n V A- AtO 

^potential ^q,~2'Y 



^potential ^2 



T 



^potential 



T 



(58) 

However, the 



We should note that the above holds even if we do not assume that j3 = a 
next equation does rely on the fact that /3 = q^. For (47a), since it is not a power function, 
we approximate the equation to first order using the Taylor expansion of e^ and show that 
the approximate equation is invariant which implies that the original equation is approximately 
invariant. This is sufficient since all of the equations in (47) are only approximations anyway. 



UJ 



UJ 




T 



frequency 



UJ 



Texp C 



^frequency 



(59) 
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Ground State Mass Profile 




Figure 6: Left: Plot of the mass profile of a ground state with its corresponding hyperbola of 
constant T overlayed. Any ground state mass profile that keeps the presented relationship with 
this hyperbola corresponds to the same value of T. Right: Examples of different ground state 
mass profiles corresponding to the same value of T. The corresponding hyperbola of constant T 
is overlayed. Notice that all three mass profiles have the same relationship with the hyperbola. 



4.2 Properties of Static State Mass Profiles 

In this last section, we discuss a few additional properties of static state mass profiles that are 
of particular interest to the study of wave dark matter. The first is the relationship between the 
total mass m and the half-mass radius for any n^^ excited state. We observe from equations 



(47c) and (47d) that the product mrh does not depend on the parameter /xq. Specifically, 



C mas sC radius fnr\\ 

mrh = ^2 ' v^^) 

where we have suppressed the notation of n. If T is constant, then, because both Cmass and 
Cradius are positive for all n (see Table [T]) , the right hand side of this equation is some positive 
constant, k, and we have that 

mrfi = k, (61) 

which defines a hyperbola. Thus, for a given n^^ excited state, all of the possible mass profiles 
for a constant value of T lie along a hyperbola. We illustrate this phenomenon in Figure [6j 

Another property of interest is the initial behavior of a static state mass profile. We explained 
previously that the mass function of any spherically symmetric solution to the Einstein-Klein- 
Gordon equations is initially cubic and approximately equal to the value in equation ([6]). It is 
routine to show that, for a static state, the value /u(t,0) = fi{0) (since the metric and stress 
energy tensor for a static state is independent of t) is 

m = M (l + ^e-2V^)^ . (62) 

Thus for an n*^ excited state and small r, we have that 

M(r> « ^ (l + ^^e--'«) . (63) 
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Ground State Mass Profile with Cubic Ground State Mass Profile with Cubic 




Figure 7: Left: Plot of the mass profile of a ground state with its corresponding initial cubic 
function overlayed. Right: Close up of the picture on the right in the region of small r. 

We illustrate this initial behavior in Figure [7j 

We make one final note here about the stability of these static states. The ground state is 
known to be stable under perturbations so long as its total mass is not too large §[5l|8}[Tll[l9] . 
In particular, the ground state is stable in the low field limit. The higher excited states do not 
appear to be stable under perturbations, instead the agitated system either settles to a ground 
state or collapses into a black hole [l|[2]. On the other hand, there has been some success in 
stabilizing higher excited states using interactions with the stable ground states |[2]. However, 
as stated in the introduction, in physical situations dark matter is always coupled with regular 
matter, which may have stabilizing effects on the dark matter. Thus understanding the stability 
of these static states in the presence of regular matter is an important open problem. It is also 
possible that dark matter may not exist as a single static state at all and as such, finding stable 
dynamical solutions of the Einstein-Klein-Gordon equations is another important open problem. 



5 Conclusion 

We summarize here the results of this paper. This paper was concerned with spherically sym- 
metric asymptotically Schwarzschild solutions of the Einstein-Klein-Gordon equations ([T]) with 
a metric of the form 

g = _e2V(t,r) ^^2 ^(^_ ?MM^ dr^ + da\ (64) 

where V and M are real-valued functions, and a scalar field of the form 

/(t,r)=e-*F(r), (65) 

where u >Q and F is complex-valued. We proved the following proposition which is designated 
here with the same number as it appears earlier in the paper. 
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Proposition 3.1 Let {N, g) be a spherically symmetric asymptotically Schwarzschild spacetime 
that satisfies the Einstein-Klein-Gordon equations ([T]) for a scalar field of the form in (14). Then 
{N,g) is static if and only if F{r) = h{r)e^"' for h real-valued and a G M constant. 

Restricting our attention, without loss of generality, to solutions where F was real-valued 
and hence the metric is static by the above proposition, we next showed that for chosen values 
of (TjHo), there was a infinite number of solutions of this type that could be distinguished by 
the number of zeros, n, the scalar field F contains. These solutions are called n^^ excited states 
when n > and ground states when n = 0. Generic examples of these static states are displayed 
in Figures [l] - |4j 

We also showed that the parameters defining such solutions, T, /ig, w, and Vq as well as the 
total mass, m, and the half-mass radius, rh, are related via the following equations 



tj"(T, Ho) ^ T exp (^Cf^equency"^^ (66a 



^0 i'^^fJ'O) ^ C potential (66b) 

m"(T,/.o)«C;;_T-3/Vy^ (66c) 
T^^iT, /.o) « a",,_T-i/Vo (66d) 

The values of the constants in these equations for various states are found in Table [T] We showed 



that these relations imply important scalings when we scale the coordinate functions as in (48) 



and restrict the types of scalings allowed if the parameter T is to be invariant (see equation 



(57])). 

Finally, we showed that the last two of the above four relations imply that for a constant 
value of T, the mass profile of any n}^ excited state lies along a hyperbola (see Figure We 
also showed that the initial behavior of a static state mass profile was cubic (see Figure]?]). 

All of these results are useful in understanding the predictions of wave dark matter in the 
case where the spacetime is static and spherically symmetric. 
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